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Abstract 

(N We have studied electron scattering by out-of-plane (flexural) phonon modes in doped suspended graphene and its effect 
b-Qon charge transport. In the free-standing case (absence of strain) the flexural branch show a quadratic dispersion relation, 
which becomes linear at long wavelength when the sample is under tension due to the rotation symmetry breaking. In 
the non-strained case, scattering by flexural phonons is the main limitation to electron mobility. This picture changes 
iy-j drastically when strains above u — 10~ 4 n(10 12 cm -2 ) are considered. Here we study in particular the case of back gate 
induced strain, and apply our theoretical findings to recent experiments in suspended graphene. 
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Graphene is a novel two dimensional material whose 
low-temperature conductivity is comparable to that of con- 
ventional metals [l(, despite much lower carrier concentra- 
tions. Interactions with the underlying substrate seem to 
be the main limitation to electron mobility, and recent 
experiments on suspended samples show a clear enhance- 
ment of mobility (more than one order of magnitude) at 
low temperatures 

In suspended graphene carbon atoms can oscillate in 
the out-of-plane direction leading to a new class of low- 
energy phonons, the flexural branch 0, 0]- I n the free 
standing case, these modes show a quadratic dispersion 
relation, so there is a high number of these low-energy 
phonons and the graphene sheet can be easily deformed 
in the out-of-plane direction. For this reason it can be 
expected that flexural phonons are the intrinsic strongly 
T-dependent scattering mechanism which ultimately lim- 
its mobility at room temperature |7fl. However, since the 
scattering process always involves two flexural phonons, 
a membrane characteristic feature, its effect could be re- 
duced, specially at low temperatures [8]. 

In the present manuscript we analyse theoretically the 
contribution of flexural modes to the resistivity in sus- 
pended graphene samples. Our results suggest, indeed, 
that flexural phonons are the main source of resistivity in 
this kind of samples. We also show how this intrinsic lim- 
itation is reduced by the effect of strain. A quantitative 
treatment of back gate induced strain where graphene is 
considered as an elastic membrane with clamped edges is 
given. 



In order to describe long-wavelength acoustic phonons 
graphene can be seen as a two dimensional membrane 
whose elastic properties are described by the free energy 

m 

J- = —k I dxdy{\7 2 h) 2 + ~ / dxdy{\u 2 vl + 2/ittf •)• (1) 



where n is the bending rigidity, A and fj, are Lame co- 
efficients, h is the displacement in the out of plane di- 
rection, and Uij = 1/2 [diUj + djUi + (dih)(djh)] is the 
strain tensor. Typical parameters for graphene [9j are 
k « 1 eV, and fi w 3 A ~ 9eVA~ 2 . The mass density 
is p — 7.6 x 10 _7 Kg/m 2 . The longitudinal and transverse 
in-planc phonons show the usual linear dispersion relation 

J^- w 2.1 x 10 4 m/s and 
1.4 x 10 4 m/s. Flexural phonons have the 



with sound velocities vl 



dispersion 
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with a — . — w 4.6 x 10 7 m 2 /s. The quadratic dispersion 



relation is strictly valid in the absence of strain. At finite 
strain the dispersion relation of flexural phonons becomes 
linear at long-wavelength due to rotation symmetry break- 
ing. Let us assume a slowly varying strain field ity(r). The 
dispersion in Eq. @ is changed to: 
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In order to keep an analytical treatment we assume uni- 
axial strain (u xx = u, and the rest of strain components 
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zero), and drop the anisotropy in Eq. ([3]) by considering 
the effective dispersion relation 



q\l a 2 q 2 + uvf 



(4) 



Long- wavelength phonons coup le to electrons in the ef- 
fective Dirac-like Hamiltonian [10| through a scalar poten- 
tial (diagonal in sublattice indices) called the deformation 
potential, which is associated to the lattice volume change 
and hence it can be written in terms of the trace of the 
strain tensor fill . Il2| 

^( r ) = 9o [u xx (r) + Uyy{v)\ (5) 

where g a « 20 - 30 eV [U . Phonons couple also to elec- 
trons through a vector potential associated to changes in 
bond length between carbon atoms, and whose compo- 
nents are related with the strain tensor as [l2 , 14 1 
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where a w 1.4 Ais the distance between nearest carbon 
atoms, /3 = 91og(t)/aiog(a) « 2 - 3 Q3, and t « 3 eV is 
the hopping between electrons in nearest carbon ir orbitals. 

Quantizing the displacements fields in terms of the 
usual bosonic a^ L ' T ' F operators for phonons of momen- 
tum q we arrive at the interaction Hamiltonian. The term 
which couples electrons and flexural phonons reads 



C-, h = EE(<+A + ) 
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where operators a t and create electrons in Bloch waves 
with momentum k in the A and B sublattices respectively. 
The matrix elements are 



K t.q,q' 
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where q = arctan (q v /q x ) and V is the volume of the sys- 
tem. The effect of screening has been taken into account 
in the matrix elements of deformation potential through a 
Thomas-Fermi -like dielectric function e (q) = 1+ e ^ Ep ' 
where V (Ep) is the density of states at Fermi energy. Note 
that g = go/sCkp) 3 eV in agreement with recent ab ini- 
tio results |13j . 

3. Resistivity in the absence of strain 

From the linearized Boltzmann equation we can cal- 



culate the resistivity as g = 



e 2 vlV(E F ) r(fe F ) ' 



where Vf 



10 6 m/s is the Fermi velocity. Our aim is to compute the 
inverse of the scattering time of quasiparticles, given by 
t^ 1 = X)k'(l — cos# k:k ')yVk,k', where Wk,k' is the scat- 
tering probability per unit time, which can be calculated 
through the Fermi's golden rule. For scattering processes 
mediated by two flexural phonons, within the quasi-elastic 
approximation, we obtain 



i=l,2 q,q' 

x n q (n q < + l)<5 k / )k -q-q'<5 (E k - E k >) (9) 



where = 1+cos 6> kk ' and / k 2 k , = 1, n q is the Bose dis- 
tribution, and i?k = vphk is t he q uasi-particle dispersion 
for the Dirac-like Hamiltonian [id ]. 

In order to obtain analytical expressions for the scat- 
tering rates it is useful to introduce the Bloch-Griineisen 
temperature Tbg- If we take into account that the rele- 
vant phonons which contribute to the resistivity are those 
of momenta q > 2kp then we have ksTsc = ftw2k F - For 
in-plane longitudinal (transverse) phonons Tbg — K 
(Tbg = 38-y/n K), where n is expressed in 10 12 cm -2 . For 
flexural phonons in the absence of strain Tbg = O.ln K. 
From the last expression it is obvious that for carrier den- 
sities of interest the experimentally relevant regime is T 3> 
Tbg, so let us concentrate on this limit. 

In the case of scattering by in-plane phonons at T > 
Tbg the scattering rate is given by [16J 
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where now g ss 3 eV is the screened deformation potential 
constant. At T <C Tbg the scattering rate behaves as 
r _1 ~ T 4 , where only the gauge potential contribution 
is taken into account since the deformation potential is 
negligible in this regime due to screening effects (r _1 ~ T 6 

In the case of flexural phonons in the non-strained case 
(in practice u -C 10~ 4 n with n in 10 12 cm -2 ), the scatter- 
ing rate at T ^> Tbg reads [l6j 
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where we have taken into account two contributions, one 
coming from the absorption or emission of two thermal 
phonons, and other involving one non-thermal phonon. 
The first one dominates over the second at T ^> Tbg- It 
is necessary to introduce an infrared cutoff frequency uj c , 
where for small but finite strain u <C 10~ 4 n(10 12 cm -2 ) is 
just the frequency below which the flexural phonon disper- 
sion becomes linear. 

From Eq. (TIT)]) we deduce a resistivity which behaves as 
g ~ T, with no dependence on n, whereas from Eq. pip 
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Figure 1: Contribution to the resistivity from flexural phonons in 
the absence of strains for two different electronic concentrations (full 
lines) and from in plane phonons (dashed line). 



we have (neglecting the logarithmic correction) g ~ T 2 /n, 
as it was deduced for classical ripples in [7|. As it can be 
seen in Fig. [31 the resistivity due to scattering by flex- 
ural phonons dominates over the in-plane contribution. 
However, this picture changes considerably if one consid- 
ers strain above 10 _4 n(10 12 cm~ 2 ), as is discussed in the 
next section. 

4. Resistivity at finite strains 

4--1- Scattering rate 

The Bloch-Griineisen temperature for flexural phonons 
at finite strains u > 1 _4 n(10 12 cm -2 ) is Tbg = 28\/Un K. 
In the relevant high-temperature regime, T 3> Tbg the 
scattering rate can be written as [I6l | 
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two terms in Eq. (|12p come from the same processes as in 
Eq. (fTTj) described above. It is possible to obtain asymp- 
totic analytical expressions for Eq. (TT2"]) . For instance, in 
the limit T < T* = ^1* ~ 7 x 10 3 u K the scattering 
rate behaves as r ~ whereas in the opposite limit it 
behaves as r -1 ~ ?L. The temperature T* characterises 
the energy scale at which the flexural phonon dispersion 
under strain Eq. ^ cross over from linear to quadratic. 

It is pertinent to compute the crossover temperature 
T** above which scattering by flexural phonons dominates 
when strain is induced. This can be inferred by compar- 
ing Eq. (JTDJ) with Eq. (JI2J and imposing tj/tf ps 1. The 
numerical solution give for the corresponding crossover 
T** « I0 6 u K. Since T** > T* we can use the respec- 
tive asymptotic expression for Eq. ([12"]) . t _1 ~ 21 to ob- 
tain T** « 32irKu/kB ~ 10 6 u K. A remarkable conclusion 
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Figure 2: Strain induced by the back gate in a suspended graphene 
membrane of length L = 1 /im as a function of the respective carrier 
density for two different AL (slack). Inset: sketch of a suspended 
graphene membrane with clamped edges. 



may then be drawn: scattering due to flexural phonons 
can be completely suppressed by applying strain as low as 
it > 0.1%. 

4-2. Back gate induced strain 

In order to compute the strain induced by the back gate 
we consider the simplest case of a suspended membrane 
with clamped edges. A side view of the system is given in 
the inset of Fig. l4~2l 

The static height profile is obtained by minimising the 
free energy, Eq. (fTJ) in the presence of the load P = e 2 n 2 /(2eo) 
due to the back gate induced electric field. The built up 
strain is related with the applied load as 0, [IH , 
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where L is the length of the trench over which graphene 
is clamped and ho is the maximum deflection (see the in- 
set of Fig. 14. 2p . We assume the length of the suspended 
graphene region in the undeformed case to be L + AL, 
where the AL can be either positive or negative. Under 
the approximation of nearly parabolic deformation (which 
can be shown to be the relevant case here [H]]) the maxi- 
mum deflection ho is given by the positive root of the cubic 
equation 



h 2 ~-LAL)h Q = 
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with trench/suspended-region length mismatch AL such 
that AL <C L. If AL = then Eq. (O can be easily 
solved and we obtain for strain 
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Figure 3: Left: Temperature dependent resistivity from Ref. [3| at 
different gate voltages; the inset shows the same in log log scale. 
Right: Result of Eq, ll6l the inset shows the back gate induced strain 
as given by Eq. II13> , 



In Fig. 14.21 the back gate induced strain is plotted as 
a function of the respective carrier density. For a typi- 
cal density n ~ 10 11 cm -2 and AL = we see that a 
back gate induces strain u ~ 10 -4 . This would imply a 
crossover from in-plane dominated resistivity g ~ T to 
g ~ T 2 /u due to flexural phonons at T** ~ 100 K, well 
within experimental reach. In the next section we will 
argue that the experimental data in Ref. [4] can be under- 
stood within this framework. Note, however, that gated 
samples can also fall in the category of non-strained sys- 
tem if AL > 0. This is clearly seen in Fig. 14.21 for AL as 
small as AL/L rj 0.3%. 

4- 3. Resistivity estimates: comparison with experiment 

Bolotin et al. |3| have recently measured the temper- 
ature dependent resistivity in doped suspended graphene. 
The experimental results are shown in the left panel of 
Fig. 14.31 in linear scale, and the inset shows the same in 
log log scale. In Ref. [4| the resistivity was interpreted as 
linearly dependent on temperature for T > 50 K. In the 
left inset of Fig. 14.31 however, it becomes apparent that 
the behaviour is closer to the T 2 dependence in the high 
temperature regime (notice the slopes of T 4 and T indi- 
cated in full lines and that of T 2 indicated as dashed lines). 
Within the present framework the obvious candidates to 
explain the quadratic temperature dependence are flexural 
phonons. Since the measured resistivity is too small to be 
due to scattering by non-strained flexural phonons we are 
left with the case of flexural phonons under strain, where 
the strain can be naturally assigned to the back gate. 

In the right panel of Fig. 14.31 we show the theoreti- 
cal T— dependence of the resistivity taking into account 
scattering by in-plane phonons and flexural phonons with 
finite strain, 



where 1/tj is given by Eq. {TUD and l/r s F tr by Eq. ljT2|) . 
We calculated the back gate induced strain via Eq. (fT5|) . 
and related the density and gate voltage as in a paral- 
lel plate capacitor model, n ~ C g {V g — V/vp)/e 0, 0| 
(C g = 60 aF//im 2 and V NP « -0.4 V). The obtained 
strain is shown in the right inset of Fig. 14.31 versus ap- 
plied gate voltage. It is seen that the system is well in 
the region where Eq. [T2] is valid. The agreement between 
left and right panels in Fig. I4.3l for realistic parameter val- 
ues [l9| is an indication that we are indeed observing the 
consequences of scattering by flexural phonons at finite, 
though very small strains. Full quantitative agreement is 
not aimed, however, since our two side clamped membrane 
is a very crude approximation to the real device 0, Q| . 

5. Conclusions 

Our theoretical results suggest that scattering by flexu- 
ral phonons constitute the main limitation to electron mo- 
bility in doped suspended graphene. This picture changes 
drastically when the sample is strained. In that case, 
strains with not too large values, as those induced by the 
back gate, can suppress significantly this source of scat- 
tering. This result opens the door to the possibility of 
modify locally the resistivity of a suspended graphene by 
strain modulation. 
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